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ABSTRACT 

We present the first results from an ongoing survey to discover and measure the metallicity of z > 3 
damped Lya systems with the Echellette Spectrograph and Imager (ESI) on the Keck II telescope. Our 
motivation arises from a recent study on the damped Lya systems suggesting only mild evolution in 
the cosmic metallicity from z ~ 2 to 4. The Echellette Spectrograph and Imager, which provides two 
complementary spectroscopic modes, is the ideal instrument for a z > 3 damped Lya survey. We describe 
our observing strategy and report on the discovery and analysis of 5 new z > 3 damped Lya systems 
acquired in a single night of observing. These observations further support the principal conclusions 
of the previous study: (1) the cosmic metallicity in neutral gas inferred from the damped Lya systems 
does not evolve significantly from z ~ 2 to 4; (2) the unweighted metallicity exhibits a statistically 
significant decrease with increasing redshift; and (3) not a single damped Lya system has a metallicity 
below [Fe/H] — —3. We discuss the implications of these results and comment on recent theoretical 
studies which attempt to explain the observations. 



Subject headings: galaxies: abundances 
1. INTRODUCTION 



galaxies: chemical evolution — quasars : absorption lines 



Determining the chemical enrichment history places fun- 
damental constraints on the processes of galaxy formation. 
Both the stellar and gas metallicity are diagnostics of a 
number of physical processes including: the star forma- 
tion rate, metal transport, the initial mass function, and 
metal production yields. For several decades, researchers 
have investigated the chemical enrichment history of the 
universe by examining in great detail the chemical evolu- 
tion of our Galaxy (Chiosi 1980; Chiappini, Matteucci, 
& Gratton 1997; Boissier & Prantzos 1999). Specif- 
ically, one examines the age-metallicity relation of vari- 
ous Galactic stellar populations and extrapolates these re- 
sults to the universe as a whole under the assumption that 
our Galaxy is not overly peculiar. This extrapolation not 
withstanding, the approach is limited by the large uncer- 
tainties in measuring stellar ages and metallicities, par- 
ticularly in the most metal-poor and presumably oldest 
stars (Edvardsson et al. 1993). While it is possible to 
obtain spectroscopy of individual stars for several nearby 
galaxies (e.g. Smecker-Hane & McWilliam 1999; Shetrone, 
Cote, & Sargent 2000), metallicities of external galaxies 
are usually determined by examining spectra of the in- 
tegrated light (Trager et al. 2000). Unfortunately this 
approach is severely limited by the well-known age, metal- 
licity, dust degeneracy and pursuing metallicity measure- 
ments beyond the local universe is extremely challenging. 
At high redshift (z - 3), Pettini et al. (2000) have in- 
troduced a method for estimating the metallicity of the 

'^Visiting Astronomer, W.M. Keck Telescope. The Keck Obser- 
vatory is a joint facility of the University of California and the Cali- 
fornia Institute of Technology. 



Lyman break galaxies, but unfortunately they have ap- 
plied it to only a single galaxy as the analysis requires 
a relatively bright galaxy in order to obtain a moderate- 
resolution, high signal-to-noise spectrum. 

Quasar absorption line studies, which examine the gas 
properties of galaxies, allow an independent and comple- 
mentary method for studying chemical evolution. In par- 
ticular, studies of the damped Lya systems - neutral hy- 
drogen gas layers with HI column densities in excess of 
2 X 10^" cm~^ - provide a measure of the mean metallic- 
ity of the universe in neutral gas from z = to 5 (Pettini 
et al. 1994, 1997; Pei & Fall 1995; Prochaska & Wolfe 
2000). Because the damped Lya systems are detected in 
absorption they are believed to present a more representa- 
tive sample of protogalaxies at high redshift (i.e. probing 
the entire galactic mass function) than galaxies observed 
in emission (Kauffmann 1996; Gardner et al. 1997; Mo, 
Mao, & White 1998, e.g.). With observations of the 
damped Lya systems, one directly studies the mean galac- 
tic enrichment as a function of time and can test scenarios 
of chemical enrichment (e.g. Pei et al. 1999). Furthermore, 
measurements of the relative abundances of elements like 
Si, Zn, Mn, Fe (Lu et al. 1996; Prochaska & Wolfe 1996, 
1997, 1999) provide an investigation into the star forma- 
tion rate and initial mass function of these protogalaxies. 

In a recent paper (Prochaska & Wolfe 2000; hereafter 
PWOO) we presented evidence that the mean metallicity of 
the damped Lya systems exhibits no significant evolution 
from z = 2 to 4. This result is in contradiction with the 
majority of chemical evolution models all of which signifi- 
cantly undcrpredict the observed enrichment at z > 3 (e.g. 
Edmunds & Phillips 1997; Pei et al. 1999; Mathhn et al. 
2000; but see Cen & Ostriker 1999; Prantzos & Boissier 
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Table 1 

NEW DLA DISCOVERED WITH ESI 



QSO 




R 


^abs 


N{m) 




[Fe/HJ 


Ion 


PSS0808+52 


4.45 


18.82 


3.114 


20.59 ±0.10 


14.15 ±0.05 


-1.94±0.11 


Fe II 2374 


PSS0957+33 


4.25 


17.59 


3.279 


20.32 ±0.08 


14.18 ±0.03 


-1.62 ±0.09 


Fe II 2344 








4.178 


20.50 ±0.10 


13.84 ±0.07 


-2.16 ±0.13 


Fe II 1608 


PSS1248+31 


4.35 


18.9 


3.696 


20.43 ±0.10 


13.89 ±0.05 


-2.04 ±0.11 


Fe II 1608 


PSS1432+39 


4.28 


18.6 


3.272 


20.95 ±0.10 


15.10 ±0.10 


-1.35 ±0.15 


Ni II 1751, Cr II 2056 



2000 as discussed below). For this epoch, however, the 
observational result was based on [Fe/H] measurements of 
only ~ 10 protogalaxies and therefore suffered from small 
number statistics. Observationally, it is very expensive to 
analyze z > 3 damped Lya systems with high resolution 
spectrographs (e.g. HIRES, UVES) because the majority 
of known z > 3 damped Lya systems lie towards faint 
quasars. To address the uncertainty arising from the small 
sample size at z > 3, wc have initiated a program using the 
new Echellette Spectrograph and Imager on the Keck 11 
telescope. The project takes advantage of the unique ca- 
pabilities of this instrument. The combination of high 
efficiency, moderately high resolution {R « 10000), and 
complete wavelength coverage from A ^ 4000 — 11000 A 
provide an ideal match to a z > 3 survey of damped Lya 
systems. In this Letter we present our first results from 
the survey, specifically metallicity measurements for five 
damped Lya systems discovered during one night at Keck 
Observatory. We compare these measurements against the 
results from PWOO and discuss the implications for chemi- 
cal enrichment in the early universe. In a future paper, we 
will present the complete set of observations and investi- 
gate chemical abundance patterns of elements like Si, Fe, 
Al, and Ni. 

2. OBSERVATIONS AND REDUCTION 

The Echellette Spectrograph and Imager (ESI; PI: Joe 
Miller, UCSC) is mounted at the cassegrain focus of the 
W.M. Keck II 10m telescope. The instrument has two 
spectrographic modes: (1) a prism-dispersed low resolu- 
tion mode with resolution R decreasing hncarly from R = 
4800 at 3900 A to i? = 1000 at 1.1/i; and (2) a multi-order 
echellette mode with R w 10000 over the entire spectrum 
(rs 11 km/s/pix). Both modes have continuous wavelength 
coverage from A = 4000 - 11000 A. 

We observed the 4 quasars listed in Table 1 on the night 
of April 7, 2000. All of them were drawn from the Palomar 
Sky Survey (Djorgovski et al. 1998) kindly made public 
by G. Djorgovski and collaborators^. At the time of the 
observations, none of the quasars were known to exhibit 
damped Lya systems. We acquired 400s exposures of each 
quasar in the low dispersion mode with a 1" slit and re- 
duced the data in real time with the IRAF task apall and 
related tasks. We then obtained two or three 1200s ex- 
posures in the echellette mode (0.75" slit lending resolu- 
tion FWHM « 45 km/s) of those quasars with promising 
damped Lya candidates; i.e., a Lya line with rest equiv- 
alent width Wx > 10 A (e.g. Wolfe et al. 1995). In this 
fashion we were able to identify a significant number of 

^http: / /www.astro.caltech.edu/ ~george/dposs/index.html 



new z > 3 damped Lya systems and measure their metal- 
licity in the same night of observing. Because the quasars 
were previously unobserved at a resolution sufficient for a 
damped Lya study, this survey will also help improve the 
statistics of the cosmological neutral gas density at z > 3 
(Storrie-Lombardi & Wolfe 2000; Peroux et al. 2000). 

We reduced the echellette observations with the echelle 
suite of IRAF tasks, in part following a recipe developed 
by S. Castro^. After subtracting the overscan and dividing 
by a normalized quartz flat, we traced the curved orders 
of a bright standard star, determined the offset for each 
faint object, and extracted the data with apall. We then 
extracted a CuArHgXe arc-lamp spectrum for each quasar 
and determined a wavelength solution using ecidentify. Fi- 
nally, the data were flux calibrated by correcting the ob- 
served flux to our observations of the spectrophotometric 
standard Hiltner 600 (Massey & Cronwall 1990), and nor- 
malized to unit flux using a package similar to the IRAF 
task continuum. In the future, we intend to reduce the 
data with a modified version of the makee packaged de- 
veloped by T. Barlow for HIRES echelle data^. For most 
of the observations, the S/N per pixel is ss 20 redward 
of Lya emission decreasing to ~ 10 at A > 9000 A and 
A < 4500 A. 

3. METALLICITY ANALYSIS AND RESULTS 

We searched each of the normalized echellette spectra 
for potential damped Lya systems and fitted the Lya 
profiles with Voigt profiles to estimate their HI column 
densities iV(HI). At a resolution R ~ 10000, the typi- 
cal Lya forest cloud is just resolved and we could account 
for the majority of contaminating features in the wings 
of the Lya profile. As an example. Figure 1 presents the 
Lya profile for the damped Lya system at z = 3.279 to- 
wards PSS0957+33. The dark curve overplotted on the 
data is a Voigt profile with logA^(HI) = 20.32 and the 
light grey lines represent a conservative la error estimate 
of ±0.08 dex. The iV(HI) values and la errors for each of 
the damped Lya systems are listed in column 5 of Table 1. 

With a FWHM resolution of « 45 km/s, we had con- 
cerns that the absorption lines from some systems would 
be unresolved and our analysis would underestimate the 
metallicity. To quantify this potential systematic error, 
we analyzed an ESI echellette spectriim of the very bright 
quasar PH957. There is a known damped Lya system at 
z = 2.3 towards PH957 with previous HIRES observations 
(Wolfe et al. 1994; Prochaska & Wolfe 1999) which re- 
veal that the metal-line profiles exhibit the majority of 

^http: / /www.astro. caltech.edu/~smc/esi. html 
* http: // spider . ipac. caltech.edu /staff/tab/ maJcee / 
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Fig. 1. — Lya profile for the damped Lyo system a,t z = 3.279 
towards PSS0957+33. The dark solid line traces the best fit value 
of Af(HI) = 20.32 and the light solid lines show an estimated Icr 
error of 0.08 dex. 



absorption within a 20 km/s interval, i.e. well within the 
ESI resolution. We compared the column density measure- 
ments determined from the apparent optical depth method 
(Savage and Sembach 1991) from 10 absorption lines ob- 
served with both ESI and HIRES. For absorption lines 
with Wx < 0.5 A all of the column density values agreed 
to within the 2a statistical error and 80% were within Icr. 
Even several obviously saturated profiles yielded column 
densities to within 0.2 dex. Therefore, we expect that 
transitions with Wx < 0.5 A will yield column density 
measurements free of saturation and we have focused on 
these transitions in our analysis. 

Depending on the absorption redshift, the HI column 
density, and the metallicity of the damped system, the 
absorber will exhibit various metal-line profiles which can 
be analyzed to assess its metallicity (e.g. Figure 2). Ide- 
ally, we measured an Fe+ column density for each system 
from the Fe II A 1608,1611,2344,2382 transitions. These 
measurements were checked for consistency against the 
7V(Si"'"), A''(Ni"'"), 7V(A1+) observations by considering the 
relative abundance patterns typically observed for the damped 
Lya systems (Lu et al. 1996; Prochaska & Wolfe 1999). 
In one case (PSS1432-h39) all of the observed Fe II pro- 
files were blended with significant absorption features or 
too weak/strong to provide a reliable column density mea- 
surement. To estimate the metallicity we relied on the ob- 
servations of other ions, in particular Ni+ and Cr+ which 
show very nearly solar relative abundances (i.e. [Ni/Fe]=0) 
in nearly every damped Lya system. For this system, we 
adopted an additional 0.05 dex error. The metallicities 
and errors are presented in Table 1 for each system. We 
also indicate the most important metal-line transition(s) 
used in determining the metallicity. 

Because we have relaxed our definition of the Fe metal- 
licity to include systems with observed Ni, Al, or Cr tran- 
sitions, we now include several other systems from our 



echelle database (BR0019-15, BRI1346-03; Prochaska & 
Wolfe 1999). Following our practice for the ESI obser- 
vations, we adopt an additional 0.05 dex error for these 
metallicities. We have also revised our [Fe/H] estimates 
for the systems towards Q1215-I-33 and QOOOO— 26 as the 
former system had an [Fe/H] value based on a saturated 
Fe II 1608 profile and we believe the Fe abundance for 
QOOOO— 26 is more reliably given by [Ni/H] which agrees 
well with the Fe abundance derived by Molaro et al. (2000). 
Finally, we also include two new observations from the lit- 
erature: the UVES observations of J0307— 4945 (Dessagues- 
Zavadsky et al. 2000) and the measurement the damped 
Lya system towards Q0201-I-11 (EUison et al. 2000). Ta- 
ble 2 summarizes all of the echelle measurements. Even 
with our relaxed criteria, we exclude a single damped Lya 
system at z = 4.203 towards BRI0951— 04 for which we 
have only detected Si II transitions, [Si/H] = —2.56 (Prochaska 
& Wolfe 1999). Because this system exhibits a low HI col- 
umn density, its inclusion would not significantly affect any 
of the following analysis except for the observed scatter in 
[Fc/H] at z > 3. 




-zoo 200 -200 200 



Relative Velocity (km s ) 



Fig. 2. — Metal-line transitions observed for the damped Lya 
system at 2 = 3.279 towards PSS0957+33. The Fe+ column density 
is well constrained by measurements of the Fe II A1608, 2344, 2832 
profiles using the apparent optical depth method. 
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Table 2 

DLA ECHELLE METALLICITIES 

QSO Zabs logiV(HI) [Fc/Hj Re? 

1 
1 
2 
1 
1 
1 
2 
1 
1 
3 
1 
3 
1 
1 
1 
4 
5 
6 
1 
4 
1 
3 
1 
4 
7 
8 
1 
6 
1 
1 
3 
6 
6 
3 
9 



"[Fe/H] estimated from Ni 
''[Fe/H] estimated from Cr 
'^[Fe/H] estimated from Al 

'^Note the tabulated metallicity for the system towards Q1223+17 is 
erroneous in PWOO. The correct value, however, was used in all of 
the figures and analysis of PWOO. 

«Key to References - 1: Prochaska & Wolfe (1999); 2: Prochaska & 
Wolfe (1997); 3: Lu et al. (1996); 4: Lu et al. (1997); 5: Prochaska 
& Wolfe (1996); 6: Prochaska & Wolfe (2000); 7: Ellison et al. 
(2000); 8: Molaro et al. (2000); 9: Dessagues-Zavadsky et al. (2000) 



4. DISCUSSION 

Figure 3 presents the; [Fc;/H] inotallicity measurements 
for z > 1.7 against absorption redshift^. The open squares 
are measurements derived from echelle observations, and 
the new ESI measurements are presented as stars. For all 
data points, the size of the data points is linear with the 

^This redshift corresponds to the optical atmospheric cutoff for 
the detection of the Lya transition 



log of the HI column density. Following PWOO, wc plot 
the column density weighted mean metallicity < Z > = 

\0g{n^etals/^gas) " log Zq = logE N{Fe+)/ ^ ^^(HI)] - 

log(Fe/H)o for two epochs: ziow = [1-5,3] and Zhigh = 
(3,4.5] with median zio^ ~ 2.20, Zhigh = 3.65. Wc have 
calculated a statistical error in < Z > based on standard 
error propagation and estimated the uncertainty due to 
sample variance with a bootstrap analysis (PWOO). We 
find < Z >iow= —1.53 ± 0.04 (statistical error) with a 
bootstrap error of 0.09 dcx and < Z >high^ —1.63 ± 0.04 
with bootstrap error 0.12 dex. Similarly the unweighted 
logarithmic means < [Fe/H] > for the two samples with 
bootstrap errors only are < [Fe/H] >iow= — 1.58±0.10 and 
< [Fc/H] >high^ -1.89 ± 0.07. Finally, the scatter in the 
two samples is (T([Fe/H]);o«; = 0.46 and a{\Pe/Yi\)high = 
0.40 each with a bootstrap error of 0.07. Comparing the 
three statistical moments for the two epochs, we note the 
values for the Zhigh sample are in excellent agreement with 
those from PWOO (there is no change for the ziow sample). 
As discussed in PWOO, there is no significant evolution 
in < Z > yet a decrease in < [Fe/H] > and cr( [Fe/H]). 
Performing a Student's t-test and the F-test on the < 
[Fc/H] > and (T([Fe/H]) statistics which test for the like- 
lihood that two the samples have a common mean and 
variance respectively, we find that they are inconsistent 
with one another at the 97% and 43% c.l. We are now 
confident that the unweighted mean < [Fe/H] > (i.e. the 
metallicity of any given damped system) is decreasing with 
increasing redshift, but less certain that the scatter in the 
[Fe/H] distribution decreases at high z. 

The additional 7 systems considered here nearly dou- 
ble the total HI content of the Zhigh sample and the new 
total {Ht = 10221 cm-2) exceeds the A^(HI) value 
measured for any single system. Therefore, the results are 
more robust to addition of individual systems. For ex- 
ample, since the publication of PWOO wc have received 
criticism that the results were largely biased by the single 

system towards PSS1443-F22 {zats = 4.23; [Fe/H] 1). 

Removing that data point from the current set of measure- 
ments reduces the mean statistics (< Z >high by 0.1 dex 
and < [Fe/H] > by 0.05) but statistically the effect is 
only significant in the scatter as this data point dom- 
inates a{\Fe/'ii\)high. Eliminating PSS1443-I-22 reduces 
a{\Fe/'ii\)high by 0.06 and a new F-test indicates the ziow 
and Zhigh distributions are inconsistent at 80% c.l. Per- 
haps more intriguing, however, is that if one were to split 
the 

Zhigh bin at z — 3.5 into two sub-samples, eliminating 
the PSS1443-I-22 system would have the effect that < Z > 
would decrease markedly (^ 0.4 dcx) at z > 3.5. We 
stress, however, that there is no physical justification for 
disregarding the system towards PSS1443+22; both the 
HI and Fe column densities are well measured and there 
may be a significant number of high metallicity systems 
at 2; > 4. As our sample size increases significantly in the 
upcoming year, we will further investigate smaller redshift 
intervals. 

Along with the results from PWOO, the new results im- 
ply an evolution in the mean metallicity of neutral gas 
which is inconsistent with the majority of chemical evo- 
lution models. This includes the analyses of Pei et al. 
(1999), Edmunds & Phillips (1997), and Malaney and 
Chaboyer (1996) all of which predict a significantly larger 



Q1331+17 


1.776 


21.18 


-2.077 ±0.041 


Q2230-F02 


1.864 


20.85 


-1.162 ±0.086 


Q2206-19 


1.920 


20.65 


-0.732 ± 0.073 


Q1215H-33'^ 


1.999 


20.95 


-1.604 ±0.122 


Q0458-02 


1.999 


21.65 


-1.642 ±0.102 


Q2231-002 


2.066 


20.56 


-1.309 ±0.100 


Q2206-19 


2.076 


20.43 


-2.634 ±0.062 


Q2359-02 


2.095 


20.70 


-1.693 ±0.103 


Q0149+33 


2.141 


20.50 


-1.799 ±0.102 


Q0528-2505'^ 


2.141 


20.70 


-1.270 ±0.140 


Q2359-02 


2.154 


20.30 


-1.905 ±0.105 


Q0216-h08 


2.239 


20.45 


-1.060 ±0.183 


Q2348-14 


2.279 


20.56 


-2.266 ±0.077 


PH957 


2.309 


21.40 


-1.804 ±0.065 


Q0841+12'^ 


2.375 


20.95 


-1.648 ±0.142 


Q2343-M2 


2.431 


20.34 


-1.228 ±0.100 


Q0201-h36 


2.463 


20.38 


-0.900 ± 0.045 


Q1223H-17^ 


2.466 


21.50 


-1.723 ±0.104 


Q0841+12 


2.476 


20.78 


-1.846 ±0.101 


Q2344-hl2 


2.538 


20.36 


-1.951 ±0.103 


Q1759-F75 


2.625 


20.80 


-1.225 ±0.108 


Q1425-F6039 


2.827 


20.30 


-1.320 ±0.058 


Q0347-38 


3.025 


20.80 


-1.828 ±0.100 


Q1055-h46 


3.317 


20.34 


-1.902 ±0.101 


Q0201-hll 


3.386 


21.26 


-1.410±0.113 


QOOOO-2619^ 


3.390 


21.41 


-2.226 ±0.142 


BR0019-15^ 


3.439 


20.92 


-1.461 ±0.159 


BRIl 108-07 


3.608 


20.50 


-2.144 ±0.101 


BRI1346-03'^ 


3.736 


20.72 


-2.663 ±0.153 


BRI0951-04 


3.857 


20.60 


-2.038 ±0.117 


BR2237-0607 


4.080 


20.52 


-2.170 ±0.167 


BRI0952-01 


4.203 


20.55 


-1.996 ±0.126 


PSS1443-I-27 


4.226 


20.80 


-0.971 ±0.114 


BR1202-07 


4.383 


20.60 


-2.220 ±0.188 


J0307-4945 


4.466 


20.67 


-1.940 ±0.333 
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Fig. 3. — [Fe/H] vs. z^bs plot of the 40 damped Lya systems which comprise our entire z > 1.7 sample. The squares indicate metaUicity 
measurements acquired from high resolution echelle observations and the stars indicate the 5 new measurments obtained with the Echellette 
Spectrograph and Imager. For two samples {zi„^ = [—1.8,3.0) and Zf^^gf^ = [3.0,4.7]) we present three statistics: (1) the Hl-weighted mean 
metaUicity < Z > represented by the solid circle with bootstrap error; (2) the unweighted mean logarithmic metaUicity < [Fe/H] > with 
bootstrap error depicted by the solid triangle; and (3) the scatter in the individual systems cr{[Fe/H]) represented by the offset vertical error 
bar. 



decrease in < Z > than observed. Prantzos & Boissier 
(2000) have noted that the absence of evolution in < Z > 
could be explained by a dust obscuration scenario where no 
system with [Zn/H] + log (HI) > 21 is observed. While 
this model matches the observations well, these authors 
did not incorporate their model within the formalism of 
Fall & Pei (1993) which is very likely to minimize the 
effects of dust obscuration. There is an ongoing survey of 
a sample of radio-selected QSO's which will help emperi- 
cally assess the effects of dust obscuration (Ellison 2000). 
Recently, IVIathlin et al. (2000) presented a comprehen- 
sive chemical evolution model for the damped Lya systems 
which included the Press & Schechter (1974) cosmological 
formalism, the Fall & Pei (1993) dust obscuration treat- 
ment, and a detailed galactic chemical evolution scenario 
along the lines of Prantzos & Boissier (2000). Their results 
have the same failing as the majority of previous models: 
an underprediction oi < Z > and < [Fe/H] > at ^ > 3. 
While they suggest that galaxy interactions and mergers 
may resolve this discrepancy, it is unclear why these pro- 
cesses would not significantly alter their results at z < 3. 
Cen & Ostriker (1999) have included the effects of merg- 
ers and feedback in their numerical simulations and their 
results indicate a milder evolution in < .Z > from z = — 2 
for objects they identify as damped Lya systems. Unfor- 
tunately, this treatment did not extend to z = 4, but pre- 
liminary results appear to agree with the results presented 



in this Letter (Cen & Ostriker 2000; priv. comm.). 

In PWOO, we noted that none of the systems with [Fe/H] > 
— 1.5 have large HI column density. This could be ex- 
plained by dust obscuration (e.g. Boisse et al. 1998), but 
it is also naturally explained by scenarios where a signifi- 
cant fraction of neutral gas has been converted into stars 
as well as galactic models with a metaUicity gradient and 
a central HI hole (Wolfe & Prochaska 1998; Efstathiou 
2000) . The latter scenarios also account for the increasing 
scatter in [Fe/H] provided the gas begins to be significantly 
consumed at z < 3. Examining Figure 3, one notes that 
none of the new systems have [Fe/H] < — 2.5 and the total 
sample [N = 19) has no system with [Fe/H] < —2.7 dex. 
We emphasize that the ESI measurements are sensitive to 
metallicities well below —3 dex^ and not a single system 
has been removed from the analysis because no metal-lines 
were detected. While future measurements may reveal a 
system or two with [Fe/H] < —3, we contend that there is 
a physical lower limit to the metaUicity of the damped Lya 
systems at all epochs. Whether this lower limit on [Fe/H] 
is a function of Pop III pre-enrichment (e.g. Wasserburg & 
Qian 2000; Ellison et al. 2000) or rapid metal enrichment 
is still a matter of debate. 

We expect our efficiency on future nights with ESI to be 

^Our 2a limits on the metaUicity are [Al/H]= -3.3, [Fe/H] = 
—2.8, and [O/H] = —4.0 assuming the feature spans 3 pixels with a 
SNR=20 and Af(HI) = 20.6 
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higher yielding at least 7 new [Fe/H] measurements at z > 
3 per night. Therefore, in the next two years our survey 
and other similar ongoing projects (Sargent, Djorgovski; 
priv. comm.) will reveal over 50 new measurements from 
ESI alone. In fact, with one more successful night we 
will have more systems in the Zhigh sample than the ziow 
sample. By pushing to lower magnitude limits {R ~ 20), 
we will also be able to examine correlations between the 
metallicity and quasar brightness. 

We acknowledge the very helpful Keck support staff for 
their efforts in performing these observations. We would 
like to thank J. O'Meara and D. Tytler for help in acquir- 
ing data on PH957. J.X.P. acknowledges support from a 
Carnegie postdoctoral fellowship and thanks ESO for their 
hospitality while this data was reduced. AMW was par- 
tially supported by NSF grant AST 0071257. 
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